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mixtures of (CE^Cl + CF^H) and (02^^012 + GPCIH2) . The 
enthalpy rond entropy values are given as isometrics from 
--40°C to 4PC)®C temperature range and I'-ISO atm. pressure 
range. For CH^Cl, G2^2^^2’ our value^ agree well with 

the older computed values hy DuPont. 



ABSTRACT 

A new vapor pressure equation lias been proposed in 
terms of polynomials in temperature wiiidli' Rave been speci- 
fically cjonstructed by a weighted Gram. -Schmidt orthiogo- 
nalization procedure. This weighting has been done with a 
weight function e“'^/^r which displays the general trend of 
a vapor pressure temperature relationship. The equation 
is as follows: 

5 = (A^ + A., X Sj, + Ag X X T?) (,) 

■ The proposed equation was tested for 50 different 
substances for which accurate vapor pressure data were 
available in the literature from T^ = 0*5 to T^ » 1, The 
results v/ere compared with the Ambrose equation (t)j the 
Martin Equation and the K.C.I). C, equation (2?)* In most 
cases Eqn. (1) was better than the Martin and K.C,D,C. 
equations. Eor Benzene, 2,2,-4trfimethyl pentane, lf20 3,nd 
CH^ our results were comparable to those of Ambrose but 
for methanol, toluene, R-butanol and oxygen Ambrose's 
equation was better. However Ambrose fitted his equation v 
only from T^ of 0,6 onwards and adjusted the number of 
constants ($n general, 8 to 10) to get the best fit. 

Also recent experimental B-Y^T and auxiliary data 
were used to compute real gas thermodynamic properties 
in the saturated and superheated regions for methyl chloride 
(OH^Cl) , difluorodichlorome thane ( CE 2 CI 2 ) » trifluorochloro-' 
methane (CE^Cl), pure compounds and for the azeotropic 
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' QMPgER 1 
INTRdDUQTrOIT 

The thesis comprises of two parts, one dealing with 

. ■¥ 

the development and testing of a vapor-pressure equation 
for organic and inorganic compounds and the other involves 
the computations of real gas thermal properties of some 
selected halomethanes and haloethanes (F,01,Br,l). 

The. development of the vapor- pres sure equation was 
prompted hy Ambrose's article(l) in -viiich he proposed an open 
ended equation for T In P in terms of independent variable 
T^ based on the Che by she v polynomials (Eg having;^ aj^ayei^ht 

function o^ 1/(1-T^)^^* lime purpose of the present investi- 
gation was to choose a wei^t function which displays the 
general behavior of vapor pressure - temperature relationship 
and using this weight function, to develop polynomials in 
temperature using the weighted Gram-Schmidt orthogonalization 
procedure (51) instead of using standard polynomials. After 
a careful study, exp(-A^/T^) was chosen as the weight fiinction. 
Hence in essence what was done was to use the following 
simple vapor- pres sure relationship^ 

Log P = B - ^ or P s= exp (B) exp(-A^/T) (l) 

and generate correction terms systematically, to lqn.( 1) 
order to represent the experimental vapor pressure - boilii^^ 
point d.ata more accurately, from the triple point to the 
critical point. The proposed equation in temas of orthogonal 
polynomials in reduced. temperature was reduced to the following 
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form ; ' : 

p' = Ca!,'"+~I2 Tj. +^31^ + Vr+— > 

Eqn.(2) truncated to 4 terms gave satisfactory results for 
approximately fifty substances tested and was found to be 
superior to tlie Ilf-K equation (27) and the Martin Equation ( 44 )* 
Ambrose equation with 8 to 10 terras was found to be better 
for some compounds than Eqn.(2). The enthalpy of vaporization 


calculated using Eqn, (2) and the Glapeyron equation 


P 


RT' 


showed good agreement with the available experimental data. 
The constants of Eqn. (2) reported in Table 2 were computed 
by the least squares orthogonal regression analysis on 
IBM 7044 Computer using selected experimejital data reported” 
in Table 1 . ■ 


In continuation with our program for the measurement 
and computation of physical and thermodynamic properties of 
Halogenated (EyCl,Br,l) methanes and ethanes (refrigerants), 
real gas thermodynamic properties such as enthalpy, entropy, etc. 
are reported here for methyl chloride (CH^Gl), diflurodichloro- 
methane (OP 2 Gl 2 )» Triflurochlorme thane (CE^Cl) and azeotropic 
mixtures (GP^Ol (50?5) + OE^H) and (0PGlH2(55.lSi) ^ - 

using the literature pressure -volme-temperature data* This 
study was primarily intended to update the existing property 
tables for GH^Gl, OE 2 CSI 2 OE^Cl using the most recent 
literature data over extended temperature and pressure rangei^ 
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Martin-Hou equation of state was used for calculating the 
departure functions (S^-S°), .etc., from 230®K to 

600 in the saturated and superheated regions. The values 
of, ^^i”^ref^’ ^®T“^ref^ were obtained using the ideal gas 
heat capacities (6^), The reference temperature was taken 
as -40°0 = -40°!’. and (S^-S^^^) were obtained by 

adding and to 

and ((Sj-S°) + (Sj - 8°^^)) respeotlfaly . 

Binary azeotropic mixture (normally minimum boiling., 
azeotrope) is found to be a 'better' refrigerant than the 
two pure components forming the mixture. The requisite data 
were available for the two mixtures (0I',C1 + CP~H) and (GJ'G1H2 
+ GI2) hence their thermodynamic properties were also 

calculated. Experimental measurement of the P-7-T data and 
other auxiliary data required for the computation of thermo- 
dynamic properties are time consuming and cumbersome. It 
would therefore be desirable to calculate the mixture proper- 
ties from the pure components using suitable mixing rules. 
Hence pure component data for GP^Cl and GE^H were taken, 
appropriate mixing rules were used (55) and vapor volumes'- and 
vapor pressures were calculated. for the ■ mixture. -These were 
C^^,iired with the expe ri|§ia.Syiial--mixtur'ej -.data, to.., check .'the.. ' . 
mlsi'i4g-c rales. If the agreement between these values is 
satisfactory, then the rmodyn^^c properties can be calculated 
using standard thermodynamic procedure. 



CHAPTER 2 


PEVBLOPMEKT QE THE .VAPOR-PRESSURE EQUATION 

The fundamental relationship relating the slope of 
the vapor pressure curve (dP/dT) and the enthalpy of vapori- 
zation (AH^) is the Clapeyron equation 


dP 

dT 


(4) 


' Tl^v ./■ W-- 

Assiiming (1) the gas is ideal (2) liquid volume very very 
small compared to the gas volume (3) represented in 

V ’ ' ' 

terms of polynomial ±n T, the following general vapor pressure 
equation is obtained. 

R 


T In P s 


A, T' 
1 




(5) 


i=i- 


The Martin equation ( '') 

/ '*• 

In P = A 4- I 4 OT + DT^ + ET^ + P log(ff-T) 

^ . ( 6 ) 

can also he represented- hy Eqn.(5) if log (S~T) is represented 

y ■ ' ' 

hy polynomial in T (1). 

Let the ' eonstants of Eqn,(5) he obtained by the least 
squares regression analysis. Mathematically the fit should 
be better when the order of the power series is. increased 
from n to n+t -vdiere n is any integer. Due to round off 
error committed during computation, _ increasing the order 
of power series from 2 to R where R iS' *t.he number of data 
points, a stage comes when increasing the order from n to 
h+1 does not improve the fit. This is because equations to 
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be solved during regression become ill-conditioned. (Tbe 

f 

equations are called ill- conditioned when they are very 
sensitive to tiny alterations caused by round off error 
in their coefficient matrix* The detailed analysis is given 
in Appendix 1). 

Recognising that pressure P can be represented in terms 
of polynomials in temperature T as in Eqn. (5), in calculating 
the coefficients of the polsmomial, the illn^conditioning of 
the coefficient matrix can be avoided through the least 
squares orthogonal polynomial fitting procedure (Appendix 1). 
According to this procedure, P is expressed in terms of 
orthogonal polynomials in as follows: 

P = ¥(Tj.) [a^ 0^ + + 

or P = W(T^).-^ . (7) 

i=o ^ ^ 

where ¥(T^) = weight function 

and 0^, 0^, 0^ 0^ are orthonormal polynomials (orthonormal 

polynomials are orthogonal polynomials with magnitude equal 
to unity) having the..follow±ag property. 


[ ¥(T^) 0^0^ dT^ = 0 


if ± ^ d 


= 1 if i =5 0 

0^, 0|, 02 “^ siTe defined as follows: 

0^ = orthogonal polynomial of 1st order = a^^ 


(9-a) 


0^ = orthogonal polynomial of Ilnd order « a^^-fa^ 2 *‘^]j. 

f9-b) 

02 = orthogonal polynomial of lllrd order a ®' 21 *'^ 22 ^r 
+^^23^r ^ (9-c)' : 


and so on 



6 


Grram- Schmidt process (Appendix 2) is used to find the 

coefficients of'Eqn,(9)» The general relationship given below 

is used to find out the orthogonal polynomials. 

n— 1 

""" “ > . (T^. <h i 

m 


T. 


0 = 
^n 




m=o 

■iir 


r^r ■" >' (T^, 0 ) 0 

m=o 


( 10 ) 


where (T^, 0 ) = inner product of and 0 = f 0 ¥(f )dT. 

* A. XHl I jC* 133. ““ 

^ (11-a) 

n-1 




m=o 


(T^, 0 )0 



i 

I f f 






m=o 


-t2 




vr(T^)dT^ 


and 


0 % 

^ 0 ,, 


W(T^) dT^ 


(11-b) 

( 12 ) 


¥2 


a 


Multiplying both sides of Sqn. (7') b 3 ’’ 0. and integrating within 

b 2 ^ ■ 

limits, all terms except ■ a. 0^ ¥(T )dT cancel because of 

. j j r r 

c3t 

the property of orthogonal polynomials -given by Eqn. (S). 
p Bij aij. = [a^ W(Tj,) alj, + ja^ 0^0^ W(Tj,)Mj, 

+ Jaj 


/ p 0 aT_ 


or 


a . 

d 


b.., 

I 0l¥(!r^) dl, 

a 


(14) 




Using Eqn. (8), equation (14.) can be reduced to the following 
expression b 




J 


P dTj, 


( 15 ) 


ill 


The values of 0^ , 0^ , 0 ^ , given^Eqn. (9) are substituted 
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in Eqn.(7) to yield 


P =s W(!r^) * ^ ^-0^00 ■*" ^^o * ^T») t ^'o * 


'11 


12 


21 


+■ + a^j Ip + - . - - 


(16) 


Eqii.(l6) can also be written by grouping the coefficients 


of as follows: 


+ Ap * + A3 * t; + A^ * TJ + . .417) 


P = ¥(T^) * |a^ 

Hence the main advantage of using the orthogonal 
polynomials is to avoid ill- conditioning of the coefficient 
matrix. A second advantage of these polsmomials is that 
the coefficients a^^a^,a 2 of the polynomials are indepen- 

dent of each other, That is the degree of approximation of 
the polynomial does not have to be chosen at the beginning 
of the computation. The values of may be computed 
successively and the decision of the number of terms to be 
used can be based on the magnitude of a^. and the overall fit 
of the data. In non-orthogonal treatment a change of the 
degree of the polynomial requires that all coefficients be 
recomputed, viidch is not desirable, 

Eqn,(l7) is the basic form of the proposed vapor pressure 
equation. The problem now is the choice of the weight function 
¥(T^) . Pitting by orthogonal polynomials primarily depends 
on'4:he selection of weight function. ¥eight function is the 
function of independent variable (temperature) and indicates 
the weightage to be given to the individual points. More- 
over it should satisfy the following criteria: 
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(a) It should he positive within the integration limits. 

(h) Its integral within the function limit should be positive. 

From Hilderbrand (3 0), the choice of the weight function 
can be made as fellows: If a function can be expressed as 
y(x) = v(x) [ 0^ 4 - 0^ + a 2 J (18) 

then the choice of ¥(x) ^-x^v(x) can be a useful one. 

One such attempt to express vapor pressure in terms of 
orthogonal pol 3 naomials was made by Ambrose who expressed 
T In P in terms of Ghebyshey polynomials. The weight function 
in case of Ohebyshev polynomials is 1/((T. -b) (a-T where 

a and b are the upper and the lower limits of the independent 
variable T,./ As weight function is the measure of weightage 
given to ind'ividual points, iimbrose (l ) selected this weight 
function with a view that as T a or T — — b, ¥(!£. ) —*- oo 
giving maximum weightage at the ends. Choice of this weight 
function gives better fit at the ends of the curve (triple 
and critical points). 

In the present study, the choice of the weight function 
was made on more realistic "a-pproach i.e. from the actual pressure 
temperature relationship. According to Ramkrishna (57) the 
best choice of the weight 'function goes to that function which 
displays the general trend of the property under investigation. 
The plot of In P vs 1/T^ is 'almost a straight line from the 
triple point to T^ = 0,8 and a S shaped curve from T^= 0,8 to 
T^ si. The linear region may be represented by 
In P - B- A^/T^ 


( 19 ) 



A 

or P = EXP(B) Bxp (- ^ ) (20) 

■^r 

Hence the weight age given to each point can be best 
represented by exp itself. Hence for the present 

study the weight function was chosen to be exp(-A^/T^). It 
satisfies the criteria number (a) and (b) as stated before, 
i.e. weight function is positive and its integration between 
integration limits to 1) is positive. 

Also Eqn. (5) can be written as 

P = Exp Bxp (Aj + ■'^3 h '^r'' ^ 

= Exp (-Ayij,) (AJ + A^ + Ajl^^ ) (22) 

Prom Hild8r]jrand' 3 criteria as given in Eqn. (18), again the 
choice of the weight function goes to expC-A^/I^) by Eqn, (22). 

Substituting W(T^) = exp(-A^/!r^) in Eqn. (17), the final 
vapor pressure equation results as given below: 

P = exp (-A^/T^) (A^ + Aj * + Aj*t2 + A^»t5+ ) 

(23) 

Eqn, (23) was fitted to the experimental data from the litera- 
ture for 50 compounds by the least squares orthogonal 
regression analysis. 
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COMPUTATIONAL PROCEDURE 


The problem Is. to compute the constants A^, a^, , 

^27 ^oo’ ^11’ ^12 Eqn. (16)* 


P = exp(-A^/T^) 


^"0 %o ■** ^"12 ^ r ^ ■** ^ 2 ^®' 21 ‘^ 22 ' 


T. 


■*^23^1*^ + ^53^r 

(16) 


3" 31 


32"r 


53"r 


3rr' 


(i) The initial value of A^ was obtained by plotting log P 
’'i^s 1/T^ and taking the slope equal to A^. 

(ii) S’qqj ®'i i» ®'12’ * • • gi''^en. by the following expressions 

(Eqns. 9, 10 and 12). 


1^0 = ® 


00 


/ 

a 


(24) 


exp (-A^/21j,)dT^ 


0. = a . . + a . _T = 
^1 11 1 2 r 


V{Po ’^r ®^PC-Ao/Tr)^''r} l*’o 

(25X 


where C = 


J— 

b 

b 


f ( 



^T^ - 


'J 

1 

J 

a 

a 


yc 

-M2 


Hence = (J0^ T^ exp(-AyT^) dT^) 0^ /0 

S'.* /n 


( 26 ) 

(27) 

"1 2 “* l/^ (28) 

Similarly for 02»03 etc,, the above integrations 

were carried out using the Romberg integration technique ( 13 ), 

In this fashion the values of a^^, ®'11» ^12* etc. were 

' ' ' . ' ■ 

calculated. Here a and b are the lower and upper limits- of 
the escperimental T^ data, 
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(iii) a^, were calculated in the following 

manner using Eqn,(l5). Hence 

j 

a 

In order to aarry out this integration, the pressure P 
was represented hy the following expression; 

P = exp (a + I + C T + P ^ B 

Using the experimental vapor pressure - boiling point data, 

constants A, B, G, etc, in Bqn,(30) were computed by the 

least squares regression analysis. Using Eqn.(30) with the 

evaluated constants, and a^^ value already computed from 

Eqn,(24-), integration of Eqn,(29) was carried out resulting 

in the value of a . 

° b 

Similarly a^ = J P(a^^ + a ^2 * ‘^^r 

Q. 

Using Eqn.(30) for P and a^^ from Eqns. (27) and (28) 

respectively, the value of was obtained. In the same 
fashion values of a^, a^ etc. were obtained. 

Thus all constants in Eqn. (16) were calculated and 
a neat expression was obtained for P in terms of 

For deciding the number of terms required to give 
the 'best* fit of the experimental data, the following 
procedure was followed. 

Eqn,(l6) was truncated after second order polynomial 
giving Bqn. (32). 

(a^^ + ^2 ^r) 
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Experimental values were substituted in Eqn*(52) and E 

calculated for each, point and compared witht he, experimental 

value. Ihe average and standard deviations of the fit were 

calculated. If the fit was not satisfactory, additional 

terms were included in Eqn,(32) and the errors were checked 

again, This process was continued upto sixth order polynomial, 

Eor most of the substances ’’best” results were obtained using 

pol 3 momials upto fourth order given below: 

a a + a . (a^ . 4 * a, « T )'+ a„(a_, + 

0 00 1 '• 1 1 12 r' ^ 2' 21 

t^23 ^r^ + ^32^r ^33 ^r ^34 ^r^j (53) 

+ ^3 ®'31 ^ 

2 


P = Exp(-A^/T^) 


^22 h 


P = Iixp(-A^/Tj,) 


(q, 3* 0,, Jt "f* ^o4 

' 0 00 1 11 ^ 2 21 


(a^a^2 + ^2^22 ^3^3:^^r ^^2®'23 ''' ^3®’33^^r ®'3®’34^ 

(34) 


Let 



a 

^00 

4* 

^1 

a.1 , 

1 1 

^2 

= ai 

^"12 

4 

^2 

CNJ 


= a^ 

^23 

4 

^3 

a M t-r 

33 

h 

= a3 

^34 





^2 ^"3 ^"31 

+ 


(35) 

(36) 

(37) 

(38) 


therefore, 

P = exp(-Ayij^) A.| + Ajlj, + Aj ll (59) 

The value of is then judiciously adjusted by increas- 
ing or decreasing the first approximation of k^ to obtain 
the ''optimum’ fit. Eqn.(39) is the proposed equation. ■ 

A program in Port ran 17 was written to compute the 
constant of Eqn.,X39) using the above procedure. The block 
diagram is given in Pigure ( 1 ), The detailed program is 
given in Appendix 3* 



BLOCK-DIAGRAM 



BIGURB 1 















0.04 


80*0 



0.04 


0.04 



Mo nof(uorobenz one 





Reduced temperature , 
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RESULTS MD DISCUSSION 


Results; Experimental data from the • literature for Oji^gen, 
water, oarbon-no no oxide, nitrous Qxlde,o*' 0 U-l'fur.. dioxide .and 
forty five organic compounds were used to test Eqn.(59)- 
The literature data references are given in Table 1, The 
critical constants were taken from Kudchadker et.al.(39). 

Eor each compound the •optimum’ values of A^ 

were calculated and these constants are reported in Table 2. 
Pointwise fit and the errors for' each compound are given 
in Table 5. The average deviation and the standard deviation 


were calculated by the following formulae: 

erroif/N 


Average deviation = 


51 . 

i =1 

Standard deviation = s;r~ 

fcr 


(error)' 


V2 


/U 


The enthalpy of vaporization, A for each compound at 
25 and at the normal boiling point (at 760 mm Hg pressure) 


were calculated using the Olapeyron equation, Bqn, (4). 

RT^ Az 




dT 


(4-2) 


The values were calculated from Eqn.(59) and 4 ^Z values 


from Haggenmacher’s equation given below 

r p,, 1 n 


z 


! 1 - 


73 I 


(43) 


The final equation used for 
Eqn^(44). 


calculations is given by 
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H = M_ _£^y2 

y 


r ( 1 - jexp(-A^/T^) (Ag/I^ + 


+ lA^/I^ 


(44) 


ODlieae calculated -iSE^ values are compared vd.th, the available 
experimental values in Table 1. 


Discussion ; In general Bqn*(39) represents the experimental 
data from the triple point to the critical point quite satis- 
factorily. further, Eqn. ( 39) was compared with the Martin ( 4 . 4 ) 
IIT-K ( 27 ), and Ambrose (I) equations. The comparison is 
reported in Table 4 in terms of the average and standard 
deviations. For most of the compounds Eqn.(39) was found to 
be superior to the Martin and IIT-K equations. Ambrose tested 
his equation ^ 

T log P = ^■"■''ag Eg (x) (45)' 

S=0 

where E (x) = Chebyshev polynomials 

X = 2T - (T ^ + T . ) / (T -T . ) 

^ max min'^ ' '• max min'^ 

JTor oxygen, nitrous oxide, methane, 2-2-4 trimethylpentane, 
benzene, toluene, methanol and n-butanol* Ambrose reported 
the deviations from the experimental value in terms of 

ln(Pgaic/^axrtasrl^ which varied from + O.OOO 4 to ± 0.0^ 

using Bqn.(45). For nitrous oxide, methane, 2-2-4 trimethyl- 
pentane, and benzene Eqn.(39) with five constants was 
comparable to Ambrose eight constant equation. However, 
for other compounds Bqn. (39) was not as good as the Ambrose 
equation. In general Ambrose used eight to ten terms in 
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Eq_ii.(45) to fit the data, Pig. 2 reports the deviations in 
Zi? = “ ^ohs compounds hy Eqn.(39), Martin 

and IIT-K equations for the region 0,8 ^ 1 .0,' Except 

for acetone, the deviations are randoni"’and smaller than 
by-.- other two models. 

Hence, a five constant vapor pressure equation has 
been proposed based on the general vapor pressure-temperature 
behaviour, which fits the experimental data from = 0,6 
upto the critical within 0.1 per cent, In general,- it is 
better than Martin's seven constant equation and III/K 
four constant equation. Ambrose equation with eight to 
ten const.ants is slightly better. Bqn, (39) can be used 
with confidence for the calculation of values and 

.for design calculations. 



CHA.PTBR 3 


REAL GAS THERMODYUAMIC PROPERTIES 


TRe literature survey on experimental pres sure -volume- 
temperature data for compounds of our interest revealed tRat 
more recent reliaRle data were availaRle for metRyl cRloride 
(01^01), difluorodicRlorometRane (CE2Ql2^* trifluorocRloro- 
metRane (CP^Cl) and for azeotropic mixtures of trifluorocRloro- 
metRane witR trifluorome tRane (CE^Cl (50^) + CP^H (50fo)) and 
monofluoromonocRlorometRane with tetrafluorodicRloroetRane 
(OPCIH 2 (55.1%) + C 2 E^Cl 2 ). The real gas tRermodynamic proper- 
ties for CH^Clj CP^Cl^ and CE,C1 Rave already been computed 
^ to 22^ ^ 

by DuPont (20) using older data. It was therefore decided to 
use the newer data and compute the real gas tRermodynamic 
properties, such as, enthalpy, entropy etc., for these compounds 
using the standard procedure (44). It is felt that these newer 
derived properties would result in "better" equipment design. 

The application of azeotropic mixture to conventional 
refrigeration systems has resulted in improved coefficient of 
performance, A knowledge of the thermodynamic properties of 
azeotropic refrigerant mixtures is thus essential to explore 
the potential of new mixed refrigerants and for the analysis 
and design of modern refrigeration systems. An azeotropic 
mixture has thermodynamic ' characteristics which for all 
practical purposes are similar to those of the pure components. 
This enables their basic properties to be computed using the 




procedure for the pure components as developed in Chapter 6. 

In ordel? to use this procedure one must carry out elahorate 
experimental investigations for the determination of P-V-T 
data and the vapor pressure data. One could avoid this hy 
developing a procedure for calculating thermodynamic properties 
of mixtures based upon the properties of pure components 
forming the mixture. One such procedure as developed by 
Piacentini and Stein (55) is tried here based upon certain 
mixing rules described in the later chapter. This procedure 
is tested for OP^Cl + CP^H mixture, 

ka a consequence of this work, a "Thermodynamic 
Package" computer program has been developed for computing 
the thermodjmamic properties of pure compounds and their mixtures, 
using the minimum computer time. 



CHAPTER 6 


GOmUTATIOIAL PROCBEURE 

The computation of real gas thermodynamic properties 
depends upon the availability of the following information: 

1 . Equation of state representing accurately, 
experimental pressure-volume-temperature data. 

2. Vapor pres sure -boiling point data 

3. liquid density vs temperature 

4. Enthalpy of vaporization as function of temperature 

5. Ideal gas heat capacity as function of temperature 

6. Eepartu re -functions calculations 

The accurac3/ of derived thermodynamic properties 
depends primarily upon the accuracy of the experimental P-V-T 
data. Por example, if the P-V-T measurements yield compressi- 
bilities accurate to 0.1 fo, then the thermal functions are 
generally accurate to about 5-10?^ because the calculation 
procedure involves differentiation twice and integration once. 
Hence one must choose an equation of state which represents 
reliable experimental data better than 0,1^ over the range of 
pressures and temperatures under investigation. 

1. After a careful and detailed study, the Martin-Hou 
equation (40) was chosen for the present study. The Martin-Hou 
equation, Eqn.(46), was derived after a systematic study of a 
large number of compounds both polar and non-polar and it has 
been successfully applied to halogenated methanes and ethanes. 
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_ 

“ v-b 


+ 


^ f^(T) fg(I) 

iW 



where f. (T) = A. + B. T + C. exp (-kT/T ) 

J. J» j_ X. o 

fg(T) = Ag + BgT + Cg exp (-kT/T^) 

f^(I) = A„ + B„T + C_ exp (-kl/T^) 


( 46 ) 


The region for each compound selected for the present 
investigation is depicted in Pig. (3). This vapor region is the 
one normally involved in the refrigeration calculations. Eqn.(46) 
truncated after five terms was found suitable to represent the 
available data with the precision and accuracy required in the 
present study. The last two teims are usually required to represen’ 
data in very high density regions. 

2 . In order to represent experimental vapor pressure- 
boiling point data accurately, the following two equations v;ere 
selected. 

log P = A + I + OT + (47) 

P = exp (-A^/T^) [ A^ + A^ T^ + A^T^ + A^T^ “ ( 39 ) 
Sqn*(47) was used because the constants A,B,C,D were already 
available in the literature. Por other compounds Eqn. (39) was 
used. 


3. Liquid densities were calculated using Yen and 
Woods procedure (^S), jt bss been demonstrated ( 68 ) that this 
procedure predicts liquid densities generally within 1 ^ and 
is convenient for machine computations. 

was calculated as follows: 


4. 
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A. + (kT +1) 0. exp(-kT ) 

(h,-h°) = Pv-RT . ( 5,3 


i=2 (i-l)(v-TD) 

(b) Entropy 


( ^.v^T ~ ^ BtK 

(58) 


(59) 


If P° is the very low pressure where gas behaves ideally, then 


v<> 


s(io^, atm) - s(P, atm) = - f (-^^)dv 

Y ^ 


( 60 ) 


R 


If gas behaves ideally also at P atm, then ( "c)P/'^ T) = . 

V\. -p V -n ^ 

Hence s° (Pjatm) - s°(p° atm) = - - dv = + - dv (61) 

0 

V V 

where s represents ideal entropy at P, atm and T.* 

Adding Eqn.(6l) to Eqn.(60)^ we get 

.0 


V 




s°(P,atm) - s(P, atm) = - ; (-^)^ dv + j | dv 

v°'- .. ^ . 

or s(P,atm) - s°(p,atm) = ^ I (-f^)^ - | I dv 
substituting Eqn.(55) into'^^Eqn. (63) we obtain v'i^oo , 

s(P,atm) - s°(P,atm) = R In (v-b) - R 3 j>(irdi 3 ) _ 


( 62 ) 

(63) 


®2 


B 


+ 


+ 


1 


3(v-b)^ 4 (Y_b)q^ ^ 2(v-b)^ 3(v-b)^ 


4 


B. 


4 - 


4 - 


C, 


4 - 


4 - 


I k 5 exp (-kT^)/T^) 


- R In V 4- 


(64) 


4(v-b) 

because v° is much larger than b and hence (v°-b)^? v°, 

(c) Real G-as Properties with respect to Reference 
Temperature of -40°C = -40°f. 

The path for calculation of enthalpy demonstrated in 
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Pig. 4 is given below; 


1. Enthalpy of saturated liquid at reference temperature 


is zero j h., = 0, 


(2) Heat the saturated liquid at ~40°C to saturated 


vapor at -40°G. 


h. = h, + Ah i 

I 1 V ! 


(65) 


' ^ '' i -4G°0 

(3) Decrease pressure of vapor from its vapor pressure 


to zero pressure at -40°0. 

h_ = - (h - h°) 


-40°C 


( 66 ) 


hp = h. + h_ = ^,- + - Ahv| -(h-h"^) (67) 

^ 11^ 1 f-40°G 

(4) Heat vapor from -40® C to temperature T at zero 


pressure 


G° d’H 
: P 


( 68 ) 


hence h^ = + h^ 2 + ^23 ~ ^1 "^^v 


-(h-ii"^) + 

1 J 


dT 


~40®C 

(69) 


(5) Increase pressure of vapor from zero pressure to 


pressure p 


h = (h^ - h°) 


(70) 


hg h^ + 11^2 + ^23 ^ ^3g 


= 


+ (hj-h°) 


■ 40®0 


^ -40°G P 


(71) 


Cj, M = logd/T^^P * 

rB,(l2-Tl^,)/2 + B5(T^-TLf.)/3 + B.(T'^-tL.)/4 (72) 


where 
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The calculation procedure of entropy at any point as 
demonstrated in Pig, 5 is given belowj 


(SjJI “ (Sji - S^) +. (s 


2^ '‘‘^2 


(s^ - 82) = #=(S2 - s'2) 


^ref ^ 


(73) 

(74) 

= (s^ -s°) + (s° - d^) r (s2-s^'*) ♦ (s°, -s°0 (75) 

Since S2, = (at low pressure) 

(Sgi-S2) = (Sj,-s^) + (s^-S2,^ (®2“®2^ 

= (s2,-s° ) + R In ^Rin -£|£ ^ (s - s°) (77) 

-= (departure function at tei'iperatare T) 

and Pressure P - R In 

— (departure function at temperature T 
and pressure P ^_) 

X c» X 

^^2”®ref^ ~ ^®2~®2^T ^^2~^ref^ 


P 

P 


,0 

"ref 


®ref 


(s2-sp^ + - - (s 


ref 


where f” 0 dT 

i P- 

T 


X log 


ref 


-s° ^ 

ref ref-^T ^ 
ref 

+ 1/T) 


ref 

(78) 

(79) 

( 80 ) 


+ B^d/T^ef . i/t2)/2 4. B^(T-T^^^) 

+ )/2 + Bg(T5-l3^p/3 


( 81 ) 

TJie Gibbs energy function can be calculated as follow 

g = h _ Is or I = I - s 

vm /— m —(s— s) 


(82) 

(83) 


\ 
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GH/kPT-ER 7 

DBSGRIPTIOIT t)F CONPUTER PROG-R/iM 

1* Main ; This program fits the P-V-T data to Eqn.(46) hy 
linear regression. Approximate values of h and k were assumed 
and average and standard errors were computed. The fit ' 
giving the least standard error was taken to be the "best" 
fit. The Main calls the subroutine Thermo. 

2, Thermo: This subroutine first fits the vapour pressure 
data. Then for any temperature and pressure, the volume of 
saturated vapor is calculated from Martin equation by Fewton- 
.■?;^p|ison technique (13). Liquid volume is calculated from 
subroutine LQDF. Heat of vaporization is calculated from 
Olapeyron equation. Since volume term is in the denominator 
inBqn. (46) it is difficult and time consuming to calculate 
V for different set of P and T values. Hence this subroutine 
reads the volimie at which v/e want to calculate the thermal 
functions. This subroutine calls EFTHA, EFTRO, HRS to. calculate 
departure fu.actio.ns, ideal gas enthalpy and ideal gas entropy 
at any temperature and pressure. Main uses Bqn, (71) and (73) 
to calculate real gas enthalpy and entropy with respect to 
datum temperature at -40°C. 

EFT HA ; This subroutine calculates enthalpy departure 
by Eqn.(57) at any pressure P and temperature T. 

4, EFTRO ; This subroutine calculates entropy departure by 
Eqn.(64) at any pressure P and temperature T. 

5, LQLF ; This subroutine calculates saturated liquid volume by 
Yen and Wood correlation (68). 


# • * • 



CHAPTER 8 


RESULTS ML LISCUSSIORS 


The sources of data used for the pure compounds and 
mixtures selected, are given in Table 5. Eor each compound, 
constants of the Martin-Hou equation, vapor pressure equation 
and the heat capacity correlation are given in Tables $> through 
The computer program listing given in Appendix lY was 
used to calculate the thermodynamic properties in the satu- 
rated and superheated regions. The same equations and the 
program were used for the two azeotropic mixtures because an 
azeotropic mixture behaves as a single component. For the 
ideal gas heat capacity of the mixture, the following formula 
was used. 


Pm 


= 0 , 




The thermodynamic properties are presented in Tables through 
for the compounds investigated at 10 deg. temperature 
interval from -40°C to + 400°C at definite selected volumes. 
Attempts were made to minimize the computational time while 
ensuring that the accuracy of results is not affected. The 
programs have been rendered flexible and are general in nature. 

The properties have been presented in units consistent with 
the international conventions. The present values are compared 
with the DuPont values in Table 1 for methylchloride , difluoro- 
dichloromethane and trifluorochloromethane. 



CHAPTER 9 


PREPICTION PROCEDURE POR MIXTURES 


Mixing Rules; 

1. Equation of States Eor Uinary mixtures, the following 
equation has been used for determining mixture constants from 
pure component constants. 

^rn ~ y-] + 2y-|y2^12 ^2 ^22 (8"^) 

'■*1 

where P^., = mixture constant 
m 

s ^2 ~ fraction of the ith component 

= pure component constants 
P^ 2 = interaction constant 

The problem is to find an expression for P^2* Piacentini 
and Stein (55 ) found that the following mixing rules for P. . 

X J 

predicted the mixture properties very satisfactorily for 
fluorocarbon system. These were adopted in the present study. 



2 . P . . 

for A^, C, 

(P^?. + - 2^.9 

' 1 - 3 . x for 


(P. . P . for 

'■11 22 


i = 3,6 
C. , i = 3,6 

i = 2,6 
k, a, T^ 


(85) 

(86) 

(87) 


for b ( 88 ) 

These rules enable the Martin-Hou constants to be determined 


= (hi + 


for a given concentration of a mixture. 
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2. Vapor Pressure Equation: The vapor pressures of 
both components of tho mixture must be fitted to identical 
model. The following mixing rule (62 ) yields results with a 
maximum error within three percent over the entire temperature 
range for some fluorocarbon mixtures: 




(89) 


Using this mixing rule, constants of Eqn.(39) for the mixture 
(OP^Gl (50?^) + CPjH) were calculated. 


3. Heat capacity of the mixture was calculated using 
Sqn. (60 ) . 

The constants of the I4artin-Hou equation, the vapor 
pressure equation, and the heat capacity equation are given 
in. Table 1 2 for CE^C1(50^) + CP^H mixture. 

Results : 

The Martin-Hou equation of state for the mixture was used 
to calculate pressure at specific temperature and volume .*and 
compared with the experimental values in Table 13 . As seen 
from this table, the agreement is satisfactory at volumes 
greater than critical volimne. Hence the selected mixing rules 
for the Martin-Hou constants are satisfactory. 

The vapor pressure equation for the mixture however, did 
not give good results. The values computed for the mixture 
are compared with the experimental values in Table 14 . It 
is therefore necessary to choose better mixing rules for the 
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vapor pressure equation constants. 

As we do not have "proper'' mixing rules for the 
prediction of mixture equations, the calculation of thermo- 
dynamic function was not attempted. 



CHA.PTER 10 


CONCLUSION 

It is hoped that (l) the proposed vapor pressure 
equation will be of use for scientists and engineers for 
the calculation of vapor pressures from the triple point 
to the critical point, (2) The calculated real gas properti 
for three pure compounds and two mixtures will prove useful 
for engineers in designing "better'’ refrigeration systems. 
The author did not have much success with the prediction 
of mixture properties from pure component data. Systematic 
work on mixing rules is necessary for this purpose. 



32 


EEPSRENGES 


1. Ambrose, D. , J.P. Counsell, and A.J. Davenport, J.Ch.ein. 

Thermo., 2 , 283 ( 1970 ), 

2. Ambrose, D. , J. Cliem. Soc . , (A) , 381 , ( 1968 ). 

3. ibnbrose , D, , B.E. Broderick, and R. Townsend, J. Chep. Soc, , 

(A), 633 (1967). 

4. iimbrose,D, , and R. Townsend, J.Chem.Soc., (A), 3614 (1963), 

5. Ambrose, D., J.D.Cox, and R.Townsend, Trans. Paraday Soc., 

56, 452 (I960). 

6 . Armstrong, G.T., P.G. Brlckwedde, and R.B. Scott, J.Res. Natl. 

Bur. Std.,^, 39 (1955). 

7. Barron, C.H., Jr., J.C. Chaty, R.A. G-onzalez, and J,¥,Bldridge , i 
J.Ohem.Eng. Data, 7, 394 ( 1962 ). 

8 . Beattie, J.A. , and S. Marple Jr., J.iAn, Chem. Soc. , 72, 1449 
(1950). 

9 .. Beattie, J.n. , S.¥. Devine, and D.R. Douslin, J.Am.Ohem, 

73, 4431 (1951). 

10. Bender, P., J.T. Purukawa, and J.R. Hyndman, Ind.Eng. Chem. , 

44» 587 (1952). 

11. Berthoud, n. , J.Chim. Phys., _1_5, 3 (1917). 

12. Brown, J.n., Ind .Eng. Ghem. , 836 (1940). 

13 . Garnahan, B. , J.n. Duther and J.O. Wilkes, Applied Numerical 
Me thods , New York, Wiley (1969). 

14. Clayton, J.O., and ¥. P. G-iauque, J.iiin. Chem. Soc . , 2610 

(1932). 

15 . Collerson, R.R., J.P. Co'unsell, R. Handley, J.J. Martin, 
and O.H.S, Sprake, J.Chem.Soc., 3697 (1965). 

16 . Connolly, J.P., J. Phys. Chem., 66 , 1082 (1962). 

17 . OouchjliJ. & K.A. Kobe,, j.Chem.Eng. Data, 6 , 229 (1961 ). 



33 


18. Crowder, G-.ii., Z.L, Saylor, T.M.Reed, and J.A. Young, 

J, Ghem. Eng. Data, _12, 481 (1964), 

19. Douslin, D.R. ,.R,S. Moore, and G. Waddington, J , Piiy , Chem . , 

63, 1959(1959), 

20. E. I. DuPont de ITemours & Company (Inc.), '* Shermo dynamic 
Properties of Metliyl Chloride” , (1939) . 

21. = E.I.Du Pont de Femours & Company (Inc.), ” She r mo dynam i c 

Properties of Difluorodichloromethane ” 1-121, (1956). 

22. S.I, Du Pont de Nemours Company (Inc,), "Ihermo dynamic 
Properties of Irifluorochloromethane ” ,1-1 3, (1959), 

23. Bucken, A., and ¥* Berger, Z. Ges. KalteADlii^, , _£l_, 145(1934). 

24. Evans, F.D, , P.E. Tiley, J.Chem. Soc., 134 (1966). 

25. Eorzia *01, A . P. , ¥. R. Norris and P.D. Rossini, J. Res. Natl. Bur. 
Std. , 43, 555 (1949). 

26. Purukawa, G.T., R.E. McCoskey, and M.L. Reilly, J. Res, Natl. 

Bur. Std., 52, 11 (1954). 

27. Goyal , V . K . , ' A. New Three Constant Vapor Pressure Equation” , 
Master's thesis, Dept ♦ of , Chem. Engg. ,IIT-Kanpur, ( 1971 ) . 

28. Gordon, J. , and ¥.P. Giauque, J.iim. Chem. Soc. , 70, 1506(1948), 

29. Hestermans, P, , and D.J. IflJhite, J.Phys . Chem. , £5, 362(1965), 

30. Hilderbrand, P.B. , Introduction to Numerical Analysis , 

Tata McGraw-Hill Publishing Co., Bombay-New Delhi, 1st 956 ) 

31. Hilderbrand, P,B. , ' Methods of Applied Mathematics ' 2nd Ed., 
Prentice-Hall (1965). 

32. Hoge, H.J., J. Res. Natl. Bur. Std., 34, 281 (1945)/ 

33. Hsu,0.C. and McKetta, J.J,, J. Chem, Eng. Data, 9( 1 ) . 45(1964). 

34. Kay, ¥.B., Ind. Eng. Chem. , 32, 358 (1940). 

35. Kay,¥. 3., J. iim. Chen. Soc, , 68, 1336 (1946). 



54 


36. Kay, ¥.3., and ¥. S. 3onliam, Chem. Eng. Sci . , 4, 1(1955). 

37. Kang, 7.L., L.J. Hirth, and J.J. McKetta, J.Chem.Eng. Data, 

6 , 220 ( 1961 ). 

38. Kobe, K.n. , H. R, Orowford, and R.V/, Stephenson, Ind.Sng. Chem. , 
47, 1762 (1955). 

39* Kudcliadker , A.?., G.H. Aland, and 3. J. Zwolinski, Ohem.Rev. 

68 , 659 (1968). 

40. Martin, J.J., Ind. Eng. Chem. , 59(12), 35 (1967). 

41. Martin, J.J,, J. Chem. Eng. Data, 5, 334 (i960). 

42. Martin, J.J, and Y.C. Hou, n.I.Ch.E.J., 5, 125 (1959). 

43. Martin, J.J. and L.3. Albright, Ind . Eng . Chem . , 188»-fi.^’ 

(1952). 

44., Martin, J.J,, "Thernodynaniic and Transport Properties of 
(J'ases, Liquids, a,nd Solids ",i-iraer. Soc . Mech.Eng., E.I., 

1959, 110 pp. 

45. McCullough, J.?., D.R. Douslin, D.¥. Scott, and R.T. Moore, 
J.ikn.Chem.Soc. , 78, 5457 (1956). 

46. McMicking, J.H., and ¥.B. Kay, 'Vapor Pressure and Saturated 
Liquid and Vapor Densities of the Isomeric Heptanes and 
Isomeric Octanes', paper presented at the 30th meeting of 
the Amer. Pet. Institute Division of Refining, Montreal, 
Canada, (1956). 

47. Michels, n. , and T‘.¥assenaar, Physica, J_4, 104 (1948). 

48. Micheli^#^^ n., I. Wassenaar., J. Walkers, C. Prins and 
L. V.D. jKlundert, J. Chem. Eng. Data, 11(4), 449 (1966). 

49. Michels, i*.. , I.Wassenaar, P. Louwerse, R.J. Lumbeck, 
and G-.I. Walkers, Physica, _19, 287 (1953). 

50. Mller, D.G-., J. Phys. Chem., ^(6), 1399 (1961 ). 

51. Murphy, K.P., , J.Chem.Eng. Data, 9, 259 (1964). 



35 


52, Oliver, G-.D, , J,¥. Grisard, and O.W, Cunningham, 
J,Am.Ohem.Soc., 22, 5719 (1951). 

53**' babotne, U. S. , H.i*, Stimson,E,I’, Hock^ and Did* bi 
J, Res. Natl, , Bur. Std. , 10, 155 (1933)* 

54. Patrick, O.R. , and G,S., Pro3ser,*-Irane,Paraday Sob*, 

60, TOO (1964). 

55. Piacentini, k» and P.P, Stein, Ohem.Bng, Prog* Symp* 
Series, 13 (81), 28, (1967). 

56. porter, P., J.Am.Ohem.Soc,, jl, 2655 (1926). 

57* Ramkrishna, D. , Ohem.Bng. Sci. (In press)* 

58i Roland.M. , Bull. Soc. Ohim., Belg* , 

59* Sage, B.H., B. C. Webster, and -W.N. ^jacey, ^ Jnd* Eng* 

: Ohem., 21, 1188 (1937). 

60* ^Selected Values of Pxroperties of Chemical Compounds” 
thermodynamics Research Center Data Project, thermo- 
: dyhai|iic Research Center, texaB' ^feMdlnivi ^ ^ College ^ 

Station, lex, (loose leaf data sheets extant | 1971). 

61. "Selected valuevO of Properties of Hydrocarbons and 
Related Compounds", iynerican Petroleum Institure Research 
Project 44, Thermodynamics Research Center, Texas A&M 
Univ, , College Station, Tex.'" ■('Loose-leaf data sheets, 
extant, 1971). 

62. Sharma,D, , "Thermodynamic Properties of a Pew Pure 
Component and Mixed Refrigerants", Master’s Thesis, 
Department of Mechanical Engineering, IIT-Kanpur.,(1969* )- 

63. Sinka, J.V., J.Chem.Eng. Data, 15(1), 71 (1970). 

64. Sinka, J.V., J. Ghem .Eng. Data, , 73 (1970). 

65. - Sokolov, V.V. , L.P. miina, and K.P. Mishchenko, Z. 

Prik. Khim. , -36 , 750 (1963). 

Suh, K.¥. and T.S. Storvick, A. I. Ch.E. J,. 13(2) . 231(1967). 


66 . 



36 


67* Tinmermans , J, "Pliysico OhenicalCo.nstants of Pure 
Organic Ooinpounds, Vol,I, *' 

Puiilisliing Company, Amsterdain, London, N.Y. ('1955). 

68. Yen, L. S.&S. S. Woods, 4.I.Cli.E., J., 12(1 ) . 95, (1966). 

69» Yesavage, V.P,, L.L. Kata« and J.E.. Powers, J.Ohein, 

Eng. Lata, H, 197 (1969). 

70. Young, S., Sci. Proc. Soc., Latlin, N.S. ZII, 374 (1910). 

71. Willingiiam, G.B. , ¥,J, Taylor, J.M. Pignocco, and 

P.L. Eossini, J. Res. Natl. Bur, Std, , 219 (1945). 

72. Zmaczymslci, M.A. , J.Chim. Phys., 27, 503 (1970), 



APPENDIX 1 


APPLICATION OP LEAST SQUARES ORTHOGONAL POLYNOMIAL PITTING 

TO ILL CONDITIONED EQUATIONS 


Tlie system of equations for which the solution is 
very sensitive to minute alterations in the coefficients 
are called ill conditioned. This tiny alteration is "brought 
ahout hy round-off error during computation, 

he the equation whose coefficients a. are to he evaluated 

2 

hy least squares regression analysis, a.’s are obtained 

J 

hy minimizing the following function f. 


N 

^ = 511 (y - ^3 

i=1 


where N = number of data points 
N 

(y - a.0,) 0. = 0 , 


df 


da. 


2 2 2 


'2 i='l 

Putting these J equations in matrix fisna , we get 
AX = B 

where A is defined as 


jJL 


'i 


A ■ = 


0101 ■ 

: 1 

0201 

1 

1 

C\J 

w 

.... ^ 


* 

> 0202 — ~ 

A. 


(A-5) 

0 -0-. 

0 1 

> 

T 






(A-2) 


(A-3) 


(A-4} 


Por 0. = x^ and for large N (Niuaher of points) and large m 
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(order of equation) and are equidistant and varying 
from 0 to 1 , then 


1 X^ 


k 

\ X dx 

j 

0 

1 

K + 1 

leads to 

r- 


- 

.w-. M 

1 

1/2 

1/3— i/ni4.1 

A = F 

1/2 

1/3 

1/4— 1/m+.2 


1/m+1 

1i/m4-2 

l/m+3 l/2m+1 


(A- 6 ) 


(A-7) 


This matrix is called Hilbert matrix and causes ill-condition- 
ing. 


The geometrical interpretation of ill- conditioning 
is given as follows; 


Let 

~ ^11 

X + a^2 y - ^1 

(A- 8 ) 


Rg = 

X + a22 y - C 2 

(A-9) 


Then the equations R^-kO and R 2=0 represent a pair of straight 
lines intersecting at point ?. Suppose, however, we do not 
demand that R^ = = 0 but only that , j R^i <C-« e 2 , 

where e^ and G 2 are small positive numbers. Then we only demand 
that the solution lie in a narrow band (bounded by the lines 
± = 0 on either side of R^ = 0 and in another narrow 

band(bounded by the lines R 2 + 62 = O) on either, side of 
R 2 = 0 , If the two lines meet at a moderately sharp angle, 
this gives us an appropriate solution. The coordinates of 
any point in the shaded parallelogram will be near to those 
of the true solution ?. 



3 ^- 


If, one the other hand, the two lines in Plg.A-l ( ) 

are nearly coincident, then the shaded parallelogram will 
he extremely long and narrow, and there will he points within 
it, situated a very great distance from P, whioii 
and | R 2 j <Ce 2 * 

4 generalization of this conce]?t provides ua with 
posaihle measure of ill-conditlohi Z£ n«3i th© eqviBttons 
represent three planes 


3 

= 0 (i*1,2,3) (A»10) 

i=i 

The smgle between any two of them, say the Ist and Ilnd is 
given hy; 

Oos, 012 ^ ^12 ^”22 ' 

' (4-11) 


Generalising this formula to n dimensions, we say that n 

linear equations represent n hyper planes and that the ( ^ 

angles between them are given by 
n 


Oos 0 




3/ . T cl . T 

ik jk 


iO 


n 


’S: 


a 


k=1 


ik 


.V2 


n 

( y- 

fer 


.2 s 1/2 
^3k^ 


(4-12) 


) 


If certain of the 0. ■ „ are very small, that is the cosines 

1 J , s 

are nearly 1, then the equations are ill-conditioned. 

Now let us see how orthogonal property is helpful in 
avoiding ill-conditioning. 
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Fron the orthogonal property. 

10^ dz = 0 if 0 7^ k 

= 1 if 3 = k 

The diagonal elements of coefficient matriz becomes 1 and 
other elements become zero. This makes the matrix unitary 
and ill- conditioning is avoided. 

Greometrical interpretation can be given as follows. 
Since planes are orthogonal, they make 90° angle between 
each other avoiding ill-conditioning as shown above by 
Eqn. (A-12). 
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APPENDIX 2 

GENEIATION OE ORTHOGONAL POLYNOMIALS BY GRAM - SCHMIDT SCHEME 

GiYen linearly independent functions 1, x, x^ 

over a certain interval |a, b| and a weight function ¥(x), 

to construct J.,0 (x)i where 0 (x) are orthogonal polynomials. 
L ) a 

^ h 


(1) 0 


1 


0 


where Ij 1 


= I ¥(x) dx 
a 


(A-13) 


!l ^il 

(2) determine C^ such that (x-C^ 0^) is orthogonal to 0^. 
(Cx - Cq0q), 0^^) = 0 (A-14) 


or 


or 


a 


(x - C^0q)0q ¥(x) dx = 0 


or X 0^^ - O^„(0, .0^) = 0 


(x,0q) 


(A-15) 

(A-16) 

X 0^ ¥(x)dx (A-17) 


0 O ' " O O ' 

(=^, 0P 

“o = 'CFTfi ■ 

X - C 0 

4 u-nwVi. 

'I 0 o;l 

In general, the process can be continued to yield, 

- S A) A 

0^ = — ^ (A-19) 


r 

X - 


> (x^,0^)0. 


m =:0 


m m 


where r takes the value of i,2,3 


r +■ 1 
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TIBLE (0 


BATA USED lOR THB a^ICULATIOST OB THERMODINAmiC PROPERTIES 


COMPOUIID ; Methyl Chloride (CH^Cl) 

( 1 ) Critical Constants ; 

T = 41 6.1 5 °E I P^ = 65.9 atra. ; y = 0.14'3t litre gia.mole'*^ 
c c c 


(2) Constants of Martin-Hou Equation of State ;( Calculated hy us) 


A^ = -0.99141179 X 10 

Bg = 0.10225509 X 

10“^ Cg =-0.2779575x1 

A^ = -0.15191194 

B^ = 0.80124886x10"^ 

C^ =0.13616686x10^ 

A^ = -0.53981769 x10"^ 

B^ = 0.11809086x10“^ 

C^=-teO. 11 11 0276x10^ 

A5 = ©.58684501x10"^ 

B^ =-0.12064974x10“^ 

O5 =0.20026005 


b = 0.004 litre mole"'’ k = 5.475 R = 0.082056 


average deviation = 0.45^ 

(5) Vaj)or Pressure Equation (35) 

log P = -6.2777215 ~ 650.82895 

T 

-0.11995288 X 10"^ X T^ 


+ 3.7672565 log T 


(4) Heat Capacity (C^) Bauation 


c 


1370.879 

T 


•1 


4.3414333 + 0.40522185 x 10 ' x T 

^5 


-0.24222874 x 10^^ x T" — -0.10767405x10^ 


Average percentage deviation = Q|3I^ 

(5) Liquid density calculated from Yen and ¥oo(b (68). 


X T' 



TABLE 7 


LATA USED FOR THE CALCULATION': OE THEBMOLYlfAMIC PROPERTIES 
COMPOUND : Difluorodichloromethane (CP2CI2) 

1. Critical Constants (39) 

T^ = 384. 7°K P„ = 38.671 atm v = 0.217 lit. g. mole"'* 

o c o 

2. Constants of Martin-Hou Equation of State (calculated By u&) 

= -0.11713538x10^ B2 = 0.15384294x10“'* 

C2 = -0.14949738x10^ A^ = 0.6770087 

B^ = -0.81274894x10“^ C^ = -0.13910429x10^ 

A^ = -0.19041529 = 0.21623801x10"^ 

C^ = 0.64101192 X 10^ A^ = 0.41118244x10"^ 

B^ = 0.12294899 xio""^ C^ =-0.55337002 x 10 

D = 0.002 lit. /mole k = +5.475 R = 0.082056 

3. Vapor Pressure Equation ; (fitted By us) 

P = Exp(-5/T^)x(-9847.6317 +32979.692 x T^ - 26721. 751x T^ 

+ 9463.4958 x t|) 

Average deviation = O.OO^'i^'p atm 
Standard deviation = 0.0001 

4. Heat Capa,city Equation : (fitted By us) 

0 = + 21.0474 + 0.1 5074x1 0"'*xT -0.12915 xIO’^^xT^ 

+ 2+ l . i |1. § - g l Q i + 0.36216 xio"^ X T^ 

Average deivation = 0,15^ 

» 

5. Liquid density calculated from Yen and Woods (68) 
correlation. 



TABLE 8 


DATi^ EOR THE GALGULiiTIOE OE THBmODYAFMIC PROPERTIES 

OOMPOBITI) ; Trif luorochlorinathane (CE^Cl) 

Critical Constants 


T^ = 379 °K = 40 atm = 0.19 lit. g. mole ' 

o c c 

Constants of Martin Hou Equation of State ; (calculated By ns) 


= -0.89171521 X 10 

B2 = 0.11677829 

xIO"”' 

02=-0.12604563x1(? 

= -0.66368465 

= 0.27641571 

xIO"^ 

C,=-0. 73918369x10 

A^^ = -0.12507003 

B^ = 0.11943893 

1 

0 

0^=0.14552154x10^ 

A^ = 0.11562561x10"^ 

3^ =-0.17963301 

x10"^ 

C^=-0. 17967315x10 


L = 0.001 lit. /mole k = 5.475 R = 0.082056 


iivarage deviation = 0.1?^ 

Va-por Pressure Equation (45) 

log^Q P(psia) = 36.7613 - - 11.80586 x log^QT 

+ 5-71495 X 10~\t 

T in °K 

Heat Capacity Equation (flttod by us) 

C = , 1 . 9 ., : . § 7 ^ 4 4 7.1 . , x , 1 of ^2,84375 - 0.2628254 x tO^^xT 

pi 

-0.35650056 x 10"^ x T^ + + 0.33124984 x 10“^xT^ 

Average deviation = 0.1?^ 

Liquid density .calculated from Yen and ¥oods(68) correlation. 



TABLE 9 


DATA USED FOR THE CALCULi;.TIOI OP THERMODYNAMC PROPERTIES 

" - n- r 'ir '.i.,' , I - r - .■in l in L. , . - -n ' ' I 'll III ' I nil — ■ , , „„ , ,,,,,, y ,,, mm ^ r™«r-rn-ir 


COMPOUND; 


Azeotropic mixture of CPCIH^ ( 55 . 1 /*^)+ 02 P^Cl 2 


Critica,! Constan cs(63) T^ = 4-15. 2°K P„ = 749 psia 

I,iiiii*.™niiii.r.rr ir rwrniiir™-. I in--Ti rr ; , ,r - 1111.111 iiin 


= 0.539 gm/cc 


Ooiistaxits of Martin-Hou Equation of State ( 63 ) 


A 2 =- 0.2351595 x 10 ^ 
= 0.3111 256 xIO^ 
A^ =- 0.1575132 xIO^ 
A5 =-0.325344 xIO"^ 

U = 0,3777405 cc/gra 


B2 = 0.2302919 xIO^ 
=- 0.2750807 x 10 ^ 
B, = 0.0 


C 2 =- 0.3169159 xIO^ 
0 ^ = 0.5292572 xIO^ 

C4 = 0,0 

0.1332969 x 10 ^ 
cc psia 


B^ = 0.201258 xIO^ 
k = 5.475 R = 0.1287070 x 10 ^ 


Vapor pressure equation ^3 ) 

log P (Psia) = 3.312656 - _ 0.771637 x 10 "^t’ 

+ 0.761974 x 10 "^xT^ 

Heat Capacity Equation, (fitted Uy us) 

Op (CPOIH 2 ) = A 2 jLii 2 | 177 i_£i 2 l + 17.978793 + 0.10870963 x 10 ”’'x T 
- 0.5347892 xio”^ xT^ + + 0.98639508 x 10 ”^xO? 


Average deviation = 0.15^ 

C^(C„P,01^) = x W ^ g 226658 x 10"^xT 

p 2 4 2 1 

-0.27010826 X 10~^ T^ + ^ + 0.94789429 x 10“^ x T^ 

T^ 

Average deviation .= 0.1^ 


Saturated Liquid Density ( 63 ) 

D(gm/cc) = 0.5391 + 0.569272 x (l-T^)^/^ + 1 . 48902 x( l-T^)^'^^ 
- 1.64833 X ( 1 -T^) + 0.988902 x ( 1 -T^)^/^ 



TABLE 1 0 


DATA USED FOR 


THE CiiLCUL/iTIOH OF THBRMODYHAMIG PROPERTIES 


COMFO JND ; Azeotropic mixt'ure of OF^Cl (505^) and 0F^H(50?^) 
Gritica.1 Constants (64 ) 

T = 292 . 6°K ; ? = 43 . 0 c atm; P = 0.564 gm/cc 

00 c 


c 

A 

A 

A 

L 


oiistants of Martin-Hou Equation of State 


2 = -9.8902899 x 10 ^ 

3 = 1.0968311 X 10 ^ 

^ = -5.8622059 x 10 
^ =-5.8740372 x 10 

= O. 344 O 8 O 9 cc/gm 


Bp = 1.2327089 
B^ =-1.0625844 
= 0.0 

B^ = 0.8136030 
k =5.475 


(64) 

O 2 = -1.8509064 
0^ = 2.9079090 

0 ^ = 0.0 
0^ = -6.4500245 
R = 0.94056886 


X lO"^ 

X 10^ 

5 

X 10'^ 

cc atm 
®K gm 


Vapor Pressure Equation (64) 

log^Q P = 8.48908 - _ 1.66749 x 10“^ x T 

+ 2.25561 X 10"^ X T^ 


Heat Capacity Equation 

0^ (CF^Cl) = z 0-3944977 . 2 , _ , x __ 10_ ^ _ 0.2628254 x 10“^ T 

p 3 T 

-0.35650056 x 10"^ x T^ + xJO^ ^ 0.33124984 x 10' 

^ . X 

0 (CF,H) = 7 1 84, 4^ , 1 IQ + 12.851214 + 0.18676057 x 10~^x T 

p .5 T 

- 0 . 1023215 X 10"^ X T^ + Ocl81926l2_x_10^ 0.21034459 x 10“®x T 




5 


Saturated liquid density (64 ) 

<3}(gm/cc) = 0.5640 + 1.2802 x (i-T^)"/^ - 1.6817 x (l-T^)^^^ 
+ 4-4210 X (1-T^) - 2.9158 x (1-T^)^'^^. 



TASLE 11 

COMPABISOK OE .Om VALUES BY-DUPOWT VALUES 


o 

o 

m 

o 




KO 


m 

VD 

cr^ 


EA 

X— 

ts^ 


"X — 

K\ 

K\ 


CO 



KA 


OJ 

04 

<M, 

CM 

CM’ 

CM 

(M 

C\J 

CM 

o 

o 


"{ — 

O 

O 

o 

O 

O 

O 

O 

o 

O 

O 

o 

o 

o 

*. 

% 

m 

<» 

■* 

« 



■* 


» 

f 

# 

o 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 


o 


o 

to 

1 — , 


m 

M3 

CM 

lA 




03 

to 



to 

T— . 


CO 

T— 

CM 

CO. 

to 

cr- 

CO 



t— 

H 

OJ 

CVl 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o 

o 

Tf— 


o 

o 

O 

O 

O 

O 

O 

O 

O 

o 

o 

o 

o 

o3 ' 

# 





* 


«■ 

« 

4e 

» 

»■ 


o 

o 

o 

o 

o 

O 

o. 

o 

o 

o 

o 

o 

o 

o 


CO 

O 


o 

o 

o 

o 

CO 

o 

o 

o 

o 

o 

cn 

o 

CO 

CA 

— 


. CO 



• 

« 


« 

« 


* 

{A 

to 


t>- 

iO 


VD 



o 

o 

o 

O 

O 

*r- 

o 

o 

CO 


VO 

CO 

LO 

<03 

CM 

1X3 

in 


• 

» 

♦ 

% 

t 

* 


CO 

CM 

CM 


in 



o 

in 

C>- 

O 

O 

CJ3 

VO 

o 

t- 

in 

CO 

T— 

CM 

VO 

o 

fen 

00 

'vT 

00 

0“ 

O 

o 

o 

trl 

r-4 

00 

o 

r— 

VD 

CM 

l>- 

O 


00 

CM 

in 

in 

MO 


* 

«• 


4 

■* 

9 ’, 

« 

% 

• 

* 

% 




C3 

in 

VD 


CO- 

in 

\o 

VD 

D- 

00 

CM 

CM 


in 

CO 

pq 






o 

o 

O 

CO 







in 

03 

00 


r- 

CM 

in 

O 

T— 

CO 

*r~ 


V- 

« 

n 

• 




♦ 

% 


« 

« 


- 

o 


00 

m 

in 

CM 

o 

CM 

CO 

T~ 

CO 

T~ 

T~ 

CM 

n 



r“ 




CM 




. 




H 

O 



O 

O 

O 

O 

O 

O 

O 

o 

o 

o 

cn 

o 

O 

> -p 

« 

€ 


♦ 

♦ 

« 

m 

• 

# 

» 

• 

• 

« 

•H 

IP 

o 

CM 

in 

in 

CM 

o 

CM 

CM 

00 


n 

CM 

o 

CM 

n 



o 

00 

o 












♦ 

t-- 

in 

CO 



o 

CO 

CM 


O- 


00 


CM 

ph a 

o 


CM 

o 


o 

EO 

G'-. 

CM 

VD 



a^ 


4J 

• 

« 

♦ 

# 


♦ 


♦ 

• 

« 


« 


« 

ci3 



VO 

c>* 



CM 

CM 

VD 

T~ 


<33 

o 

n 





T~- 



t — 


CM 

Y— 



CM 

T“ 






CM 















1 — 1 










H 





O 






O 




r3 

o 

o 

o 

o 

CM 

o 

o 

O 

O 

o . 

EO 

o 

O 

o 

bd 

o 

o 

tn 


tn 

in 

O 

o 

in 


l>- 

O 

o 

EH o 

EO 



O 

CM 

to 

EO 



O 

CM 

to 



o 


CM 

CT\ 


O 

m 



